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Abstract.  The  BIPM/IAU  Joint  Committee  on  relativity  for  space- 
time  reference  systems  and  metrology  (JCR)  has  worked  in  collaboration 
with  the  IAU  Working  Group  on  relativity  for  celestial  mechanics  and  as¬ 
trometry  (RCMA)  on  the  problems  of  astronomical  relativistic  space-time 
reference  frames  and  the  RCMA  has  proposed  to  the  IAU  a  consistent 
framework  for  defining  the  barycentric  and  geocentric  celestial  reference 
systems  at  the  first  post-Newtonian  level.  Because  new  clock  technology 
and  space  missions  will  necessitate  the  application  of  this  framework  for 
time  and  frequency  measurements  in  the  solar  system,  practical  resolu¬ 
tions  are  proposed  by  the  JCR  to  the  IAU  to  allow  a  consistent  treatment 
of  these  measurements.  The  paper  outlines  the  conclusions  of  the  work 
which  resulted  in  the  proposed  IAU  resolutions. 


1.  Introduction 

In  its  1991  General  Assembly  the  International  Astronomical  Union  (IAU)  ex¬ 
plicitly  adopted  the  general  theory  of  relativity  as  the  theoretical  framework 
for  the  definition  and  realization  of  space-time  reference  frames  (IAU,  1991). 
Barycentric  and  geocentric  coordinate  time  scales  and  the  relativistic  transfor¬ 
mations  between  them  were  defined,  together  with  procedures  for  their  realiza¬ 
tion. 

Previous  work  in  the  CCDS  Working  Group  on  the  application  of  gen¬ 
eral  relativity  to  metrology  (Guinot,  1997  and  references  therein)  has  explicitly 
applied  the  IAU  1991  framework  to  time  metrology.  When  studying  the  conse¬ 
quences  of  the  increasing  accuracy  in  the  realization  of  the  SI  units,  and  particu¬ 
larly  of  the  SI  second,  it  recognized  that  the  realized  and  expected  progresses  in 
time  metrology  were  such  as  to  require  an  extension  of  the  IAU  1991  framework. 
Based  on  these  conclusions  and  recognizing  that,  in  this  context,  metrology  could 
not  be  separated  from  the  definition  and  realization  of  space-time  reference  sys¬ 
tems,  the  BIPM  proposed  to  the  IAU  to  create  a  Joint  Committee  to  study  these 
matters. 

At  the  Kyoto  General  Assembly  (1997),  the  IAU  endorsed,  by  its  Resolution 
B3  (1997),  the  creation  of  the  Joint  Committee  on  Relativity  for  space-time  ref¬ 
erence  systems  and  metrology  (JCR).  Also  the  IAU  Working  Group  on  Relativity 
in  Celestial  Mechanics  and  Astrometry  (RCMA)  continued  under  the  chairman¬ 
ship  of  M.  Soffel  and  a  common  document  was  established:  a  “List  of  issues” 
(document  available  on  the  web:  www.bipm.fr/WG/CCTF/JCR/jcrissue.html) 
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discussing  as  much  as  possible  all  topics  that  may  be  addressed  by  the  two  work¬ 
ing  bodies.  Circulars  describing  the  work  of  the  JCR  are  issued  as  required  and 
the  web  pages  of  the  JCR  (www.bipm.fr/WG/CCTF/JCR)  are  updated  with  each 
Circular. 

At  an  early  stage,  some  discussions  took  place  on  whether  the  formalism 
should  be  based  strictly  on  the  General  Relativity  theory  or  on  some  parametrized 
theory.  Although  no  complete  consensus  could  be  reached  at  that  time,  the  de¬ 
tailed  work  started  with  General  Relativity  on  the  basis  that  it  was  more  likely 
to  reach  a  final  stage,  i.e.  a  form  that  can  be  adopted  as  a  Recommendation  to 
extend  the  IAU  1991  framework.  The  RCMA  initiated  a  theoretical  document 
aiming  at  describing  completely  the  reference  systems,  the  coordinate  transfor¬ 
mations,  and  the  equations  of  motion  in  the  solar  system.  The  JCR  initiated 
studies  on  the  shortcomings  in  the  definition  and  usage  of  the  IAU  1991  frame¬ 
work  as  well  as  on  the  application  of  the  extended  framework,  with  particular 
emphasis  on  time  and  frequency,  time  transformations  and  the  realization  of 
coordinate  times. 

In  section  2,  we  recall  the  content  of  the  IAU  1991  resolution  A4  dealing 
with  General  Relativity  and  expose  some  of  the  limitations  of  this  framework. 
The  practical  case  of  the  extension  of  the  IAU  1991  framework  for  time  and 
frequency  applications  in  the  solar  system  is  discussed  in  section  3,  which  forms 
the  basis  for  the  new  resolutions  to  be  proposed  to  the  IAU  and  other  concerned 
bodies. 


2.  The  IAU  1991  framework 

2.1.  The  IAU  1991  recommendations 

The  IAU  resolution  A4(1991)  contains  nine  recommendations,  the  first  five  of 


which  are  directly  relevant  to  our  discussion. 

In  the  first  recommendation,  the  metric  tensor  for  space-time  coordinate 
systems  (/,  x)  centered  at  the  barycenter  of  an  ensemble  of  masses  is  recom¬ 
mended  in  the  form 

9oo  — 

_1  +  2f/,,2'x)+o(o-n 

9ot  = 

CS 

1 

CO 

(1) 

9ij  = 

6,,  (i  +  AI/T)  +  o(c-). 

where  c  is  the  speed  of  light  in  vacuum  (c  =  299792458  m/s),  U  is  the  sum  of 
the  gravitational  potentials  of  the  ensemble  of  masses  and  of  a  tidal  potential 
generated  by  bodies  external  to  the  ensemble,  the  latter  potential  vanishing  at 
the  barycenter.  The  algebraic  sign  of  U  is  taken  to  be  positive.  This  recommen¬ 
dation  recognizes  that  space-time  cannot  be  described  by  a  single  coordinate 
system.  The  recommended  form  of  the  metric  tensor  can  be  used  not  only  to 
describe  the  barycentric  reference  system  of  the  whole  solar  system  resulting  in 
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the  BRS,  but  also  to  define  the  geocentric  reference  system  (GRS)  centered  in 
the  center  of  mass  of  the  Earth  with  a  suitable  function  U,  now  depending  upon 
geocentric  coordinates.  In  analogy  to  the  GRS  a  corresponding  reference  system 
can  be  constructed  for  any  other  body  of  the  solar  system. 

In  the  second  recommendation,  the  origin  and  orientation  of  the  space  co¬ 
ordinate  grids  for  the  solar  system  (BRS)  and  for  the  Earth  (GRS)  are  defined. 
Notably  it  is  specified  that  the  space  coordinate  grids  of  these  systems  should 
show  no  global  rotation  with  respect  to  a  set  of  distant  extragalactic  objects.  It 
also  specifies  that  the  SI  (International  System  of  Units)  second  and  the  SI  meter 
should  be  the  physical  units  of  proper  time  and  proper  length  in  all  coordinate 
systems.  It  states  in  addition  that  the  time  coordinates  should  be  derived  from 
an  Earth  atomic  time  scale. 

The  third  recommendation  defines  TCB  (Bary centric  Coordinate  Time) 
and  TCG  (Geocentric  Coordinate  Time)  —  the  time  coordinates  of  the  BRS 
and  GRS,  respectively.  The  recommendation  also  defines  the  origin  of  the  times 
scales  (their  reading  on  1977  January  1,  0h0m0s  TAI  ( JD  =  2443144.5  TAT) 
must  be  1977  January  1,  0h0m  32.184s)  and  declares  that  the  units  of  measure¬ 
ments  of  the  coordinate  times  of  all  reference  systems  must  coincide  with  the  SI 
second  and  SI  meter.  The  relationship  between  TCB  and  TCG  is  given  by  a 
full  4-dimensional  transformation 


TCB  -  TCG  =  c 


-2 


-y  +  Uext(t,xE(t)) 


dt 


vErE 


0(c 


-4\ 


(2) 


where  x%E  and  vE  are  the  barycentric  coordinate  position  and  velocity  of  the 
geocenter,  rE  =  xt  —  xlE  with  xl  the  barycentric  position  of  the  observer,  and 
Uext(t,XE(t))  is  the  Newtonian  potential  of  all  solar  system  bodies  apart  from 
the  Earth  evaluated  at  the  geocenter. 

In  the  fourth  recommendation  another  time  coordinate,  Terrestrial  Time 
(TT),  is  defined  for  the  GRS.  It  differs  from  TCG  by  a  constant  rate  only 


TCG  -TT  =  Lgx  (JD  -  2443144.5)  X  86400,  LG  «  6.969291  X  10“10,  (3) 

so  that  the  unit  of  measurement  of  TT  agrees  with  the  SI  second  on  the  geoid. 
TT  represents  an  ideal  form  of  TAI ,  the  divergence  between  them  being  a  con¬ 
sequence  of  the  physical  defects  of  atomic  clocks.  The  fifth  recommendation 
states  that  the  former  barycentric  time  TDB  may  still  be  used  where  disconti¬ 
nuity  with  previous  work  is  deemed  to  be  undesirable. 

2.2.  Some  problems  in  the  framework  of  the  IAU  1991  recommen¬ 
dations 

Because  of  the  form  of  the  metric  (1)  in  the  IAU  1991  framework,  time  trans¬ 
formations  and  the  realization  of  coordinate  times  in  the  barycentric  system  are 
not  specified  at  the  c-4  level,  i.e.  at  a  level  of  a  few  parts  in  1016  in  rate.  Within 
that  approximation,  constants  Lg  and  Lq  were  introduced  in  notes  to  the  IAU 
1991  Recommendation  3  to  express  the  mean  rates  between  time  scales  as 
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TCB  -  TDB  =  Lb  x  {JD  -  2443144.5)  X  86400, 

Lb  ~  1.550505  X  10“8,  and 

TCB-TCG  =  Lc  X  (  JD  -  2443144.5)  X  86400  +  vlErlE/c2  +  P, 

Lc  ~  1.480813  X  10“8,  (4) 

where  P  represents  periodic  terms.  Since  JD  is  not  specified  to  be  a  particular 
time  scale,  these  constants  were  not  properly  defined  so  that  confusion  appeared 
in  their  usage.  Following  Fukushima  (1995)  we  propose  to  use  as  more  rigorous 
definitions,  in  order  to  overcome  the  above  problems  <  TCG/TCB  >=  1  —  Lc 
and  <  TT/TCB  >=  1  —  Lb,  where  <>  means  a  sufficiently  long  term  average 
taken  at  the  geocenter.  The  actual  computation  of  Lc  and  Lb  requires  the 
integration  of  solar  system  ephemerides  and  the  specification  of  an  averaging 
duration.  This  process  may  be  applied  to  the  utmost  accuracy,  after  a  choice  of 
ephemerides  and  averaging  duration.  For  example  Irwin  and  Fukushima  (1999) 
determined  Lc  =  1.48082686741  X  10-8  ±  2.  X  10-17.  However,  because  no 
unambiguous  definition  may  be  provided  for  Lb  and  Lc,  these  constants  should 
not  be  used  in  formulating  time  transformations  when  it  would  require  knowing 
their  value  with  an  uncertainty  of  order  1  X  10-16  or  less. 

Another  problem  arising  from  the  situation  before  IAU  1991  is  that  the 
barycentric  dynamical  time  TDB  did  not  have  a  good  definition.  This  could 
be  corrected  by  turning  a  specific  value  of  Lb  into  a  defining  constant  thus 
providing  in  retrospect  a  good  definition  of  TDB.  It  was  not  felt  necessary  to 
address  explicitly  this  point  in  a  recommendation. 

The  most  important  questions,  however,  concern  the  extension  of  the  frame¬ 
work,  at  an  uncertainty  level  beyond  that  covered  by  the  1991  formulation.  This 
means  to  apply  the  extended  post-Newtonian  metric  (S  off  el,  2000)  to  the  realiza¬ 
tion  of  the  barycentric  coordinate  time  TCB  and  to  the  time  transformations  at 
an  adequate  level  of  uncertainty,  and  is  described  in  section  3.1.  Finally,  some 
consequences  of  the  application  of  the  IAU  1991  framework  in  the  geocentric 
system  will  be  considered  in  section  3.2. 

3.  Time  and  frequency  applications  in  the  solar  system:  basis  for 
proposed  Resolutions 

For  practical  applications  concerning  time  and  frequency  measurements  in  the 
solar  system,  it  is  necessary  to  consider  a  conventional  model  for  the  realization 
of  time  coordinates  and  time  transformations.  This  model  should  be  chosen  so 
that  i)  its  accuracy  is  significantly  better  than  the  expected  performance  of  clocks 
and  time  transfer  techniques,  ii)  it  is  consistent  with  the  general  framework  of 
(S  off  el,  2000)  and  iii)  it  may  readily  be  used  with  existing  astrometric  quantities, 
e.g.  solar  systems  ephemerides. 

Regarding  item  (i),  we  may  derive  reasonable  accuracy  limits  in  a  straight¬ 
forward  way.  At  present  the  best  accuracies  are  reached  by  Cs-fountain  clocks 
at  less  than  two  parts  in  1015  in  fractional  frequency  (Simon  et  al. ,  1998), 
with  a  frequency  stability  characterized  by  a  standard  Allan  deviation  of  or¬ 
der  cry{T)  =  4  X  10-14r-1/2,  for  an  integration  time  r  in  seconds.  In  the  near 
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future,  high  accuracy  laser  cooled  Rb  clocks  (Bize  et  al. ,  1999)  and  space-borne 
Cs  clocks  (Salomon  et  al.,  1999)  are  expected  to  reach  accuracies  of  a  few  parts  in 
1017  in  fractional  frequency  and  stabilities  of  order  CFy{r)  =  1  X  10  14r  1/2 ,  with 
t  in  seconds.  The  uncertainty  in  the  time  transformations  should  induce  errors 
that  are  always  lower  than  the  expected  performance  of  these  future  clocks.  In¬ 
cluding  a  factor  2  as  safety  margin,  we  therefore  conclude  that  time  coordinates 
and  time  transformations  should  be  realized  with  an  uncertainty  not  larger  than 
5  X  10-18  in  rate  or,  for  quasi-periodic  terms,  not  larger  than  5  X  10-18  in  rate 
amplitude  and  0.2  ps  in  phase  amplitude. 

For  the  spatial  domain  of  validity  of  the  transformations,  we  note  that 
projects  like  SORT  (Solar  Orbit  Relativity  Test)  plan  to  fly  high-accuracy  clocks 
to  within  0.25  AU  of  the  Sun,  but  it  turns  out  that  formulations  valid  there  will 
also  be  valid  at  a  few  solar  radii  from  the  Sun  which  is  therefore  the  lower  limit 
for  the  distance  to  the  barycenter  that  we  will  consider.  In  the  geocentric  system 
we  will  consider  locations  from  the  Earth’s  surface  up  to  geostationary  orbits 
(|X|  <  50000  km). 

To  comply  with  item  (ii),  we  explicitly  derive  the  developments  in  the  next 
sub-section  from  the  general  framework  presented  in  the  RCMA  report  (S  off  el, 
2000),  using  the  same  conventions  and  notations.  Regarding  item  (iii)  we  show 
below  how  the  time  transformations,  e.g.,  TCB  — TCG ,  may  be  performed  with 
the  existing  astrometric  quantities  and  tools. 

3.1.  Barycentric  celestial  reference  system 

Following  (S  off  el  2000),  the  metric  tensor  in  the  BCRS  (this  is  the  new  name  of 
the  BRS,  where  C  stands  for  Celestial)  is  expressed  as 

goo  =  -(1  -  K^o^x)  +  wL{t,x))  +  |-(wo(l,x)  +  A(/,x))), 

got  =  -^wt{t,x),  (5) 


where  (t  =  TCB,x)  are  the  barycentric  coordinates,  wq  =  G^aMa/ta,  with 
the  summation  carried  out  over  all  solar  system  bodies  A,  r^  =  x  —  xa,  ta  =  | r ^4 1 , 
and  where  wl  contains  the  expansion  in  terms  of  multipole  moments,  as  defined 
in  (S  off  el,  2000),  required  for  each  body.  The  values  of  masses  and  multipole 
moments  to  be  used  may  be  found  in  IAU  or  IERS  documents  (IERS  1996),  but 
care  must  be  taken  that  the  values  are  in  SI  units  (not  in  so-called  TDB  units 
or  TT  units).  The  vector  potential  wl(t,x)  =  ^raj(fx)  and  the  function 
A(/,x)  =  J2a  ^a(Tx)  are  given  later. 

From  (5)  the  transformation  between  proper  time  and  TCB  may  be  derived. 
It  reads: 


1 


1 


dr/dTCB  =  1 — 5-  w0  +  wL  +  —  +—  [ --v4  - -v2w0  +  4vlwl  + -w^  +  A 


1 


1 


(6) 
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To  reach  the  uncertainty  mentioned  earlier,  it  is  sufficient  to  express  the 
vector  potential  wA(t,x)  of  body  A  as 


wA(t,  x)  =  G 


~(r4  X  S^)8 

2ri 


Mav\ 

rA 


(7) 


where  is  the  total  angular  momentum  and  vA  is  the  coordinate  velocity  of 
body  A.  As  for  the  function  AA(t,x)  it  is  sufficient  to  express  it  as 


A  i(f,x) 


GMa 

rA 


2  Gv\ 


rkAakA 


(8) 


where  rEA  =  |X_B  —  xa|  and  aA  is  the  coordinate  acceleration  of  body  A. 

Evaluating  the  A^-terms  for  all  bodies  of  the  solar  system,  we  find  that 
|A^(t, x)|/c4  may  reach  at  most  a  few  parts  in  1017  in  the  vicinity  of  Jupiter 
and  about  1  X  10-17  close  to  the  Earth.  Presently,  however,  for  all  planets  except 
the  Earth,  the  magnitude  of  A^(t,x)/c4  in  the  vicinity  of  the  planet  is  smaller 
than  the  uncertainty  in  wq/c 2  or  wE/c2  originating  from  the  uncertainties  in  its 
multipole  moments  so  that  it  is  practically  not  needed  to  account  for  these  terms. 
Nevertheless,  when  new  astrometric  observations  allow  us  to  derive  the  moments 
with  adequate  uncertainty,  it  will  be  necessary  to  do  so.  In  any  case,  for  the 
vicinity  of  a  given  body  A,  only  the  effect  of  A A(t,x)  is  needed  in  practice,  i.e., 
the  effect  of  ^b^Ab(1,x)  is  smaller  than  our  accuracy  specifications.  For  a 
clock  in  the  vicinity  of  the  Earth,  to  be  compared  with  other  clocks  in  the  solar 
system  or  to  TCB ,  it  may  thus  be  needed  to  account  for  Ae(f,x)/c4.  Also,  in 
these  formulas,  the  terms  in  are  needed  to  reach  the  mentioned  uncertainty 
only  for  Jupiter  ( S  ~  6.9  X  10 38m2s~1kg)  and  Saturn  (S  ss  1.4  X  10 38m2s~1kg), 
in  the  immediate  vicinity  of  these  planets. 

Similarly,  the  transformation  between  TCB  and  TCG  may  be  written  as 


TCB 


TCG  =  c“2 


^ v2Ew0ext(yLE ) 


v2  \ 

+  w0ext(yLE)  j  dt  + 


4vEwlxt(xE)  +  -wlext{yiE, 


dt 


—  (3wQext(xE)  +  n|;/2)t>^r^], 


(9) 


where  t  is  TCB  and  where  the  index  ext  refers  to  all  bodies  except  the  Earth. 
This  equation  is  composed  of  terms  evaluated  at  the  geocenter  (the  two  integrals) 
and  of  position  dependent  terms  in  rE,  with  position  dependent  terms  in  higher 
powers  of  rE  having  been  found  to  be  negligible.  The  first  integral  may  be 
computed  from  existing  planetary  ephemerides  (Fukushima,  1995;  Irwin  and 
Fukushima,  1999).  Since,  in  general,  the  planetary  ephemerides  are  expressed 
in  terms  of  a  time  argument  Teph  which  is  close  to  TDB ,  rather  than  in  terms 
of  TCB ,  the  first  integral  will  be  computed  as 
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It  (if  +  w°ext^EIj  dt 


[  (~7r  +  W0 eii(xfi)]  dteph 

'teph0  \  1  J 


/(I- Lb).  (10) 


Terms  in  the  second  integral  of  (9)  are  secular  and  quasi-periodic.  They 
amount  to  ~  1.1  X  10-16  in  rate  ( dTCB/dTCG )  and  primarily  a  yearly  term  of 
~  30  ps  in  amplitude  (he.  corresponding  to  periodic  rate  variations  of  amplitude 
~  6  X  10-18).  Terms  in  A ext(t,  x)  that  would  appear  in  this  integral  are  negligible 
in  the  vicinity  of  the  Earth  where  this  formula  is  to  be  used.  Position  dependent 
terms  in  c-4  (the  last  two  terms  in  (9))  are  not  negligible  and  reach,  for  example, 
an  amplitude  of  0.4  ps  (~  3  X  10-17  in  rate)  in  geostationary  orbit. 

The  material  in  this  sub-section  is  expressed  in  the  Resolution  ‘Extended 
relativistic  framework  for  time  transformations  and  realization  of  coordinate 
times  in  the  solar  system’  proposed  to  the  IAU  General  Assembly  in  August 
2000. 


3.2.  Geocentric  reference  system 

Evaluating  the  contributions  of  the  higher  order  terms  in  (Soffel,  2000),  it  is 
found  that  the  IAU  1991  framework  with  the  metric  of  the  form  (1)  is  sufficient 
for  time  and  frequency  applications  in  the  GRS  in  the  light  of  present  and 
foreseeable  future  clock  accuracies.  In  particular  a  choice  of  gauge  and  of  the 
state  of  rotation  of  the  GRS  (kinematically  or  dynamically  non-rotating)  does 
not  influence  the  time  and  frequency  application  at  the  given  level  of  accuracy. 
However,  some  care  needs  to  be  taken  when  evaluating  the  Earth’s  potential  at 
the  location  of  the  clock  which  is  not  trivial  when  accuracy  of  order  10-18  is 
required  (Klioner,  1992;  Wolf  and  Petit,  1995;  Petit  and  Wolf,  1996). 

Presently,  the  time  scale  of  reference  for  all  practical  matters  on  Earth  is 
Terrestrial  Time  TT  or  one  of  the  scales  realizing  it  and  differing  by  some  time 
offset  ( e.g .,  T AI,UTC,GPS- time).  The  definition  of  TT  reads:  “time  scale 
differing  from  the  Geocentric  Coordinate  Time  TCG  by  a  constant  rate,  the 
unit  of  measurement  of  TT  being  chosen  so  that  it  agrees  with  the  SI  second 
on  the  geoid.”  According  to  the  transformation  between  proper  and  coordinate 
time,  this  constant  rate  is  given  by  d(TT)/d(TCG)  =  1  —  Ug/c2  =  1  —  Lq,  where 
Ug  is  the  gravity  (gravitational  +  rotational)  potential  on  the  geoid. 

The  gravity  potential  on  the  geoid  Ug  is  determined  with  an  uncertainty 
slightly  lower  than  1  m2  s-2  (Groten,  1999)  so  that  the  corresponding  rate  un¬ 
certainty  in  the  definition  of  TT  is  of  order  1  X  10-17.  In  view  of  the  perfor¬ 
mance  of  future  clocks,  it  may  be  desirable  to  change  the  definition  of  TT  in  a 
way  that  would  ensure  continuity  but  would  not  be  limited  by  this  uncertainty. 
One  solution  is  to  define  TT  differing  from  TCG  by  a  fixed  constant  rate,  i.e. 
Lq  becomes  a  defining  constant.  To  ensure  the  best  continuity,  this  conven¬ 
tional  value  of  Lq  should  be  chosen  as  the  current  best  estimation  of  Ug/c 2,  i.e. 
Lq  =  6.969290134  X  10-10,  from  the  value  Ug  =  62636856  m2  s-2  provided  by 
the  International  Association  of  Geodesy  Special  Commission  3  (Groten,  1999). 

The  material  in  this  sub-section  is  expressed  in  the  Resolution  ‘Re- definition 
of  Terrestrial  Time  TT’  proposed  to  the  IAU  General  Assembly  in  August  2000. 
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